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It has been suggested that major transitions in evolution require the emergence of novelties, often 
associated to the cooperative behaviour of previously existing objects or agents. A key innova¬ 
tion involves the first cooperative interactions among molecules in a prebiotic biosphere. One of 
the simplest scenarios includes two molecular species capable of helping each other forming a cat¬ 
alytic loop or hypercycle. The second order kinetics of the hypercycle implies a hyperbolic growth 
dynamics, capable of overcoming some selection barriers associated to non-cooperative molecular 
systems. Moreover, it has been suggested that molecular replicators might have benefited from a 
limited diffusion associated to their attachment to surfaces: evolution and escape from extinction 
might have been tied to living on a surface. In this paper we propose a field theoretical model of 
the hypercycle involving reaction and diffusion through the use of a many-body Hamiltonian. This 
treatment allows a characterisation of the spatially correlated dynamics of the system, where the 
critical dimension is found to be dc = 2. We discuss the role of surface dynamics as a selective 
advantage for the system’s survival. 

PACS numbers: 05.40.-a, 89.75.-k, 64.60.-i, 87.15.Zg 


I. INTRODUCTION 

Life on our planet has experienced several key innova¬ 
tion events since its early appearance [T] . The emergence 
of self-replicating molecular entities pervaded the rise of 
complex life forms, including those having an embodied, 
protocellular organisation. Each of these so called ma¬ 
jor transitions incorporated some kind of novel way of 
enhancing cooperation among simpler subsystems. Most 
transitions are essentially tied to improving the selective 
advantage of the agents inhabiting the new complexity 
level, and that means exploiting some class of coopera¬ 
tive interaction at the higher scale. In particular, it has 
been conjectured that, in order to achieve higher levels 
of genetic information, moving beyond small, unreliable 
chains of molecular units, cooperation among different 
replicating molecules was required m- 

One of the most celebrated candidates to explain this 
class of phenomena and its implications is Eigen’s theory 
of the hypercycle [H [S] that proposes a key role of mu¬ 
tually enhancing replication among complementary sub¬ 
units. Specifically, the general hypercycle involves a pop¬ 
ulation of n coupled molecular species described by a set 
of concentrations {pi} with i = l,...,n whose replica¬ 
tion is dependent upon the presence of another species 
that act as catalyst, forming a closed loop. Typically, 
the deterministic dynamics is given by a set of Hrst-order 
coupled equations 


^=A,(p;{r,,}) (1) 
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with p = (pi,...,p„) and where {Ti^} is a set of pa¬ 
rameters. In its original formulation, these equations 
included a reaction and an outflow term, i. e. = 
Uli{p) — piF{p). Here Uli{p) includes the reactions among 
molecular species and the last term is a continuous dilu¬ 
tion flux F required to maintain a constant population 
c = Pj of molecules. Using c = 1, it is not difficult to 
show that the dynamics is described by a general form 


dfh 

dt 


i=i 


( 2 ) 


Which automatically satisfies the so called constant pop¬ 
ulation constraint (CPC), i. e. dc/dt = 0. 

The specific form of the original n-member hypercy¬ 
cle is a closed loop of catalytic reactions. This can be 
described by the set of equations: 


dpi p 

^ i-ji—lPiPi—l Pi / ^ ^ j,j—^PjPj — ^ 


( 3 ) 


where the i-th member of the cycle catalyses the growth 
of the {i l)-th species while it is also helped by the 
i — 1-th one in the loop (with the constraint j -|- 1 —>■ n 
when i = n ) [MI]- 

This model has several remarkable properties. One is 
that it properly describes the logic of molecular cooper¬ 
ation networks m where cross-catalytic molecules are 
present (particularly RNA or peptide chains). Secondly, 
the hypercycle involves a superexponential kinetics that, 
under ideal circumstances predicts a finite-time singu¬ 
larity. An important implication of this growth dynam¬ 
ics is that the hypercycle is capable of overcoming ex¬ 
ponentially growing, Darwinian replicators [2]. Finally, 
although we will restrict our discussion to simple, molec¬ 
ular replicators, it is worth mentioning that it has also 
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been used as a relevant descriptor of other out of equilib¬ 
rium systems, including ecosystems mm and economic 
networks m- Similarly, it has been suggested that in a 
homogeneous (well mixed) setting this cooperative struc¬ 
ture would easily break down (and become extinct) under 
the presence of parasites [15] . 

The emergence of the first hypercycles is a specially 
relevant problem. If we consider the problem of repli¬ 
cating molecules experiencing mutation and selection, a 
key result of theoretical models of evolution is that there 
are sharp limits to the complexity that such system can 
achieve m- For a given mutation rate, it can be shown 
that there is a limit to genome size (scaling up as the 
inverse of mutation rate). How can this be overcome? 
One way is to make several genomes cooperate, thus ef¬ 
fectively moving beyond the so called error threshold. In 
an early biosphere, before closed compartments might 
have facilitated chemical reactions, pairwise interactions 
among the members of the hypercycle must have been 
difhcult. Stochastic effects and molecular decay might 
counterbalance the potential for rapid growth predicted 
by the ideal model. Not surprisingly, several studies have 
focused on the role played by space in the evolution and 
persistence of these structures (Mil]. 

The simplest model considers two classes of replicating 
chains m which can self-replicate provided that a pair¬ 
wise interaction takes place (see Fig0. This requirement 
implies that a cross-catalytic reaction is needed and a 
second-order kinetics underlies the population dynamics 
of both chains. Such kind of reaction-diffusion dynamical 
system can be naturally explored using field theoretical 
approximations, which allow to understand the univer¬ 
sal features of certain classes of reaction kinetics and the 
role played by spatial degrees of freedom. This approach 
has been successful in approaching a wide range of prob¬ 
lems including chemical kinetics, population dynamics or 
molecular processes both within cells and in cell-free con¬ 
texts m as well as to the dynamics of quasispecies HHj. 
This is also the approach we take here to investigate some 
universal properties of hypercycles. 

The paper is organised as follows. In section II we 
outline the basic results associated to the mean field de¬ 
scription of a two-member hypercycle. In section III a 
master equation is introduced on a discrete spatial lat¬ 
tice and make use of Doi-Peliti’s formalism to obtain a 
hamiltonian description of the system. Then the contin¬ 
uum limit is taken such that a held theoretical action may 
be written. We perform a perturbative one-loop compu¬ 
tation around a non-trivial vacuum corresponding to the 
stable active state. We will hnally apply renormalization 
group ideas in order to discuss the couplings flow equa¬ 
tions on a coarse-graining scheme. This leads to a unique 
insight on the role of dimensionality for the spatial dy¬ 
namics of the system. In particular, we will be able to 
observe these effects by plotting a phase diagram of ac¬ 
tive/inactive states on the large scale or infra-red limit. 
In section IV we briefly discuss potential extensions to 
Wdimensional hypercycles. In section V we summarise 



FIG. 1: Bifurcations in the simple symmetric hypercycle. The 
logic organisation of this pairwise cooperative loop is outlined 
in (a). In (b) the fixed points of the system are plotted against 
the control parameter F (here K = 1 and S = 0.125). A first 
order transition is at work. For F < Fc an absorbing state 
(extinction of both cooperators) is present, whereas a stable 
hyper cycle can form if F > Fc and the initial state is such 
that p(0) > P-. 


our basic results and discuss possible further work. 


II. MEAN FIELD THEORY OF THE 
SYMMETRIC HYPERCYCLE 

The simplest model of a hypercycle involves a two- 
member system with symmetric interactions, as sum¬ 
marised in Fig 1. In this section we make the assumption 
of a well-mixed scenario, where the population densities 
of each member of the cooperative loop, hereafter indi¬ 
cated as A and B, will be given by pi and p 2 , respec¬ 
tively. As shown in the diagram, each molecular species 
replicates with help from the other partner while both 
degrade. 

The hypercycle equations for this system would read, 
in general, as a pair of coupled cooperation-like differen¬ 
tial equations: 

f (4, 

^ = r2lPlP2 (^1 - - (52P2 • (5) 

The coefhcients F^j stand for the replication rate of each 
molecular species under the presence of the second one. 
There is a limit to the maximum population of molecular 
replicators provided by the carrying capacity K. The 
last term stands for linear degradation rates. It is worth 
mentioning that this model includes several particular 
cases, some of which have been studied in previous papers 
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FIG. 2: The potential function associated to the mean field 
theory for the symmetric hypercycle, as given by ([^. Here we 
fix if = 1 and 5 = 0.125 which gives a critical point at Fc = 1 
(as in figure 1). The parameter F takes different values: (a) 
F = 0.5, (b) r = Fc = 1.0, (c) F = 1.25 and (d) F = 1.25, 
respectively. 


An alternative approach to stability is obtained by as¬ 
suming that the dynamics is derived from a potential 
function, which allows us to write our system as 


dp _ dV{p) 
dt dp 

which in our case reads 





( 9 ) 


( 10 ) 


This function is such that its minima and maxima cor¬ 
respond to the stable and unstable fixed points, respec¬ 
tively. Four snapshots of this landscape are displayed in 
Figj^ It can be noticed that the single-well potential is 
deformed as we approach F^ and two minima appear after 
we cross it, with the valley associated to the stable hy¬ 
percycle becoming deeper for larger F, while po becomes 
a metastable state of the system [2T1 [22]. 


III. FIELD THEORY 


flO] . For example, if if 1 and r 2 i = 0, 1^2 = 0 we 
recover a model by Schnerb et al Bl¬ 
under the symmetry assumption, ri 2 = r 2 i = F and 
Si = $2 = S, equal populations are achieved at equilib¬ 
rium, so the fixed points are, Pi = P 2 = P- A. good 
description of the dynamics [21 is provided by the single 
differential equation model for the density p: 

( 0 , 

The fixed points of this equation are the trivial one, po = 
0, associated to the extinction (absorbing) phase plus two 
solutions 



where we use Fc = 86/K. The (linear) stability of these 
fixed points is determined by the sign of 


^r(Ps) 



( 8 ) 


for each ps. Stable and unstable points are characterised 
by negative and positive values of A, respectively. 

The marginal stability condition Ar = 0 provides the 
location of a first-order phase transition, as displayed in 
Figj^ The trivial fixed point is stable under small pertur¬ 
bations, whereas p_ and p_|_ are unstable and stable, re¬ 
spectively. But these two nontrivial points only exist for 
F > Fc, and they collide at Fc so that pi(rc) = p(rc) = 
if/4. Here, a sharp transition governed by a saddle-node 
bifurcation takes place. As expected, if degradation is 
negligible {5 —)■ 0) then p+(0) —)■ if/2. 


In this section we approach the problem of stochastic 
hypercycles along the lines of previous work involving a 
field theory analysis of nonequilibrium systems, particu¬ 
larly for reaction-diffusion chemical systems [2HIl3()j . In 
order to address the stochastic, spatial behavior of our 
two-member hypercycle, we consider a spatially-extended 
system that is discretized as a lattice in d dimensions. 
Our goal is to define the conditions under which the hy¬ 
percycle is expected to be stable and avoid the absorbing 
state and determine the role of dimensionality. 


A. Reaction-diffusion model 

Let us consider an L'^ lattice of cell size 1. In each lat¬ 
tice site we allow the system to evolve under the second- 
order, catalytic reactions 

A-f H-^4 (2A-k B, A-f 2H) (II) 

as well as the (linear) decay transitions 

(A,B)4^0. (12) 

Here, p stands for the reaction rate of the catalytic pro¬ 
cesses of A via presence of B and viceversa. We also let 
both species undergo diffusion at equal rate D. In addi¬ 
tion to the reactions above, one needs to introduce extra 
processes that render the hard-core repulsive forces [28j . 

(2A + B, A+ 2B) ^ (2A, 2B) (13) 

and, similarly, we have 

(2A-kB, A + 2B) — 


A + B. 


(14) 
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This can be shown to be equivalent to van Wijland’s con¬ 
struction [25] for hard-core diffusive particles. Later, we 
will set A = A' for simplicity. Now the potential func¬ 
tion becomes bounded from below and it is possible to 


recover the first order phase transition predicted in the 
previous section (see Fig§. Spatial dynamics can be 
implemented via a master equation for all lattice sites 
r,sG (L/Z)'^} 


d 

—P{fi,fh;t) = Djius + 1)P(..., nr - 1, Us -I- 1,..., m; t) - nrPjn, rh; t)] 4- 

{r,s} 

S[{nr + 1)T’(..., rir -I- 1,..., to; t) — nrP{n, m; t)] + 

r 

y^ n[{nr — l)mrP {..., rir ~ 1) ■ • ■) to; t) — nrmrP{n, to; t)] -|- (15) 

r 

y^ X[nr{nr — l){mr + l)P{n ,..., rrij. -|- 1,...; t) — nr{nr — l)mrP{n, to; t)] + 

r 

y^ X'[{nr + l)nr(rnr + 1)T’(. ■ ■ ,rir + 1, ■ ■ ■, rrir -|- 1,...; t) — nr{nr — l)mj.P{n, to; t)] -|- {n GG to , } 

r 


where D is the diffusion coefhcient, r is summed over 
all the lattice sites, and s over the nearest neigbors of 
r. The first line in the RHS of (15) implements diffusion 
through a random hopping of particles between nearest 
neighbor sites. The rest of the lines generate the inter¬ 
actions from (11)-(14|. Initial conditions, P{n,rfi,0), are 
typically chosen as a Poisson distribution, with an aver¬ 
age density per site equal for both types of particles. 


B. Doi-Peliti second quantisation approach 


Following the Doi-Peliti second quantisation procedure 
[?7lE5] . we construct a field theory for two species of 
diffusive molecules A and B and write down an Euclidean 
action upon which we will be able to apply perturbation 
theory. 

We introduce two sets of creation and anihilation op¬ 
erators at each lattice site. 


ttr \nr} = n \nr — 1) bj. |toj.) = TO |to^ — 1) 
al \nr) = \nr — 1 ) bl \mr) = \mr — 1 ) 

which fulfill the standard commutation rules 


L, bl 


= Sr 


(16) 


By incorporating these rules, the operators have a 
bosonic character. Using the vacuum state |0), defined 

by 


Ur |0) = hr |0) = 0 

we construct an orthonormal basis of states |n, to) for the 
Fock space, defined by 


\n,m)=l[{alr’-{blr^\0), (17) 

r 

In terms of this Fock space, the state of the system is 
given by the vector state \P{t)), defined as 


\Pit)) ='^P{n,mA)\n,rn). (18) 

n,rn 

In terms of this vector state, the master equation is analo¬ 
gous to a Schrodinger equation in imaginary time, namely 

||P(t)) = -i7|P(t)), (19) 

with a Hamiltonian H = Hu -\- Hs + + H\ + Hy 

defined by diffusion terms plus interactive terms 

at)(a,. - a,) -k (bl - bl)(br - S^)! 

{rs) 

Hs = |(aj - l)ar + (bl - l)6r| 

r 

Hf_i = ^'^^(al+ bl — 2)alblarbr . ( 20 ) 


For the sake of simplicity, we skip Hx and Hx' terms. 
Eq. (19) can be formally solved in terms of the operator 
H yielding 


\P(t)) = eM-Ht) |P(0)). 


( 21 ) 


This defines a Schrodinger equation that is real and 
thus the approach departs from the standard many-body 
quantum mechanics. In particular, the states are linear 
functions of the probabilities, instead of the amplitudes 
of probability. 
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FIG. 3: Computation shown in (a) is associated to the propagator renormalization. Which couples to the mass (decay) coupling 
renormalization. On the other hand, (b) is the r 2,2 vertex computation, that allows the one-loop renormalization of the coupling 
/i. Time flows to the left. Simple lines stand for the gaussian propagator of a-fields, Gq and double lines for 6-fields, Gq. 


C. The Hypercycle action 

From the previous discrete set up, it is posible to realize 
a continuum limit considering that the scale at which we 
are interested to perform measures is much larger than 
the size of the lattice cells. In fact, this is equivalent to 
introducing a physical cutoff A ^ y. 

In order to achieve a field theoretical description, a 
coherent state representation is introduced by setting 


recovered when one-point expected values are taken and 
fluctuations are ignored. Thus, we read off the associated 
vacuum expected values (vev) for each field: {(f)) = (^/>) = 
1, while, at the active phase, {(j)) = ('0) = p+. 

Since we are interested in probing the system fluc¬ 
tuations around the active stable solution (see Figj^), 
all fields need to be shifted to their corresponding min¬ 
imum before implementing weak coupling perturbative 
RG methods: 


dr l^r) — 4‘r \4^r} \'^r} — '^r I'^r) 

{(j)r\al = {(j)r\ (j)* = {tprl'tp* ■ 

This is chosen such that the scalar product is 
((/)* |(/)r) = , this structure will allow a resolution 

of unity so that it will be possible to treat (p* and ip* 
as auxiliary fields independent from (f>r and ipr, respec¬ 
tively. Thus, after taking the continuum limit. 


(p —y u -[-1 
Ip ^b+1 


(p —^ a -l- p-^- 
tp ^ b-\- p+ . 


This shift leads us to the a modified field theory action 
]. Now, the interaction terms 


(221 with a, a, b, b fields 
obtained are 


{(pr} -t (p{x, t) {tpr} -t '(p{x, t) 

{ipl} ^ , 

the system is mapped into a statistical field theory m- 
Hence, the action describing the two-member reaction- 
diffusion hypercycle is given by 


Uo 

Hi 

n2 


—Di^aV'^a + W^6) -I- S(aa + bb) + Sp(a + b) 
—p{a + b}(a + 1)(5 -|- l)(a -b p)(b + p) 


b{a 1)^(0 + p) + o,(b l)^(fe + p) (fl + p){b T p) 


H. = X' 


a{a + p) + b(b -bp) (a -b 1)(6 -b l)(a -b p)(6 -b p). 


Slp^cp^ip^tp] 


r r ~ ~ ^ 

/ d‘^x / dt^pdtp + tpdtpJ + ( 22 ) 

i=0 


Here onwards p+ is replaced by p for simplicity. Com¬ 
puting the classical solution by functional analysis on the 
system’s action, i.e., 


The interaction parts "Hi, i G {0,..., 3} correspond to 
the reactions (11|-(14), respectively. The MF solution is 





6=0 


= 0 
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and ignoring all fluctuations (n-point correlation func¬ 
tions), we recover the MFT equations 0 and (0 with 
X = X' = K-\ 


D. Perturbative expansion and RG flow 


div? 

db 


= 2m 
£9 + 


2 _ ^ 

IT 1 + m 
2 9^ 

TT 1 + 


2 ’ 


(23) 

(24) 


Firstly, let us suppose that the steady microscopic con¬ 
figuration is such that the system sits somewhere over the 
tipping point of Fig. Thus, we would like to compute 
how fluctuations can drive the system down to its inac¬ 
tive absorbing phase and figure out under which condi¬ 
tions the stationary state of the system will remain active 
as we let the scale flow towards the infra-red (IR) limit. 

Via straightforward power counting we can identify 
that there is a dimension below which "Hi operators be¬ 
come marginally relevant. Introducing a diffusive tem¬ 
poral scaling [t] = where k is the momentum scale, 
and taking D as an adimensional constant we have 

[a] = [b] = [p\ = , [d] = [b] = 


Henceforth, for d = 2 the coupling constant p, becomes 
marginal. Note, that the system is not posed near a 
critical state, thus, we cannot extract universal proper¬ 
ties from a perturbative field theoretical computation. 
Nonetheless, it is possible to implement a coarse-graining 
process at the vecinity oi d = 2, i.e., e = 2 — d and work 
up to one-loop expansion using the Feynman diagrams 
shown in Figj^ This scheme will ultimately lead to the 
construction of a phase diagram (Figj^). 

Let us note that for the shifted action ( [2^ a new mass 
term is obtained, namely, := S — pp + 2>Xp^ [55] . 
This mass coupling plays the role of an effective decay 
rate for the active system and it lays on the denominator 
of the fields propagators. From the shifted action (22) 
we derive an effective action valid only in the vicinity 
of d = 2 by only inculding 'Ho and Hi interactions. A 
total of 26 vertices are obtained. After introducing A as 
the large momentum cutoff, i.e., A ^ 1//, we define the 
adimensional couplings 


= A-^M^/D , 
n = A~‘^p , 
9 = A-^pjD . 


The analysis of these flow equations allows us to probe 
the dynamics of the system at different length scales. 
Lines in Figj^are solutions to these equations given dif¬ 
ferent initial conditions, which are understood as the mi¬ 
croscopical coupling constant values. As b grows the cou¬ 
plings are readjusted to a new value that renders the dy¬ 
namics of the new and larger length scale. From these 
results we read how the system’s activity pervades under 
rescaling. 

The sign of the e term in (24) shifts the nature of the 
trivial fix point from unstable in d < 2 (Figj^) to d = 3 
(Figj^) for the coupling g. 

Solutions plotted in Figj^are restricted from a pertur¬ 
bative point of view [51], thus, the active/inactive phase 
diagram provided by Figj^ is not exact once a certain 
threshold is exceeded. Nonetheless, the criterion used to 
select the curves that drive the system into extinction is 
a conservative one. If the mass coupling flows above the 
values shown in Fig. 4a and 4b, then we assume that 
the system falls into extinction, eventhough at the strict 
& —>■ oo limit that curve might return to a zero mass 
value. 

In this approach, a radical transition in the shape of 
the active/inactive curve shown in Figj^ is promptly ob¬ 
served when allowing the dimensionality go from 2D to 
3D, suggesting a clear advantage for the system’s sur¬ 
vival on a 2-dimensional layout. This observation is a 
direct result of the field-theoretical formalism. This is 
the main result of our paper, indicating that the survival 
of a cooperative molecular system (as the one described 
here) is strongly enhanced under a dynamics constrained 
on a surface. This can be easily interpreted in terms of 
the well known conjecture that reaction kinetics might 
have been strongly facilitated in two-dimensional sub¬ 
strates where limited movement of molecules is involved. 
This would not be the case in a three-dimensional con¬ 
text, where a high effective probability of extinction is 
expected. 


IV. N-MEMBER HOMOGENEOUS 
HYPERGYCLE 


The one-loop perturbative computation is derived from 
the vertex-functions Fi.i and r 2,2 which diagramatic ex¬ 
pansions are shown in Fig. These calculations generate 
coarse-grained couplings ga and m|j. Upon differentia¬ 
tion with respect to the momentum scale k the couplings 
flows are obtained. 

Next, we introduce a length-scale parameter b = log(s) 
with s = I'/I, where I' is the coarse-grained scale. Thus, 
6 ~ and 6 —> oo will correspond to the IR limit. The 

results are given by: 


In this section we will briefly outline how dimension¬ 
ality affects a generalization of hypercycles known as the 
N-member perfectly cyclic systems |5] using the same 
field-theoretical setup as in the previous analysis. We 
have already defined in (3) the basic deterministic equa¬ 
tions for such a general hypercycle. Extensive work has 
shown that complex dynamics can emerge, including os¬ 
cillations and chaos |S1|55HS7|- Is it possible to generalise 
our previous approach to an arbitrary hypercycle? In 
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FIG. 4: Coarse-graining scale flow towards the IR {b —> oo) 
of the couplings for the dimensionality d < 2 (a) and d = 3 
(b). The shaded areas correspond to an inactive phase on the 
large-scale limit. Figure (c) corresponds to a superposition 
of (a) and (b). Here a drastic shift is observed as we let the 
system undergo rescaling from a 2D setup to a 3D one. In 
particular, mortality becomes predominant as d is increased. 
The numerical analysis has been carried out using an RK4 
algorithm. For numerical simplicity, on the y-axis g' = {2/n)g 
is plotted instead of g. Since n > 1/4, figures (a) and (b) are 
obtained for n = 1/4. 



FIG. 5: Higher-order hypercycles can be build by the cooper¬ 
ation of multiple replicators forming a closed loop. Here we 
show a homogeneous N = 6 hypercycle, where each member 
helps the next one. Decay transitions are not displayed. 


this section we aim to address this problem using the 
field theory formalism. 

Consider a perfectly ordered symmetric hypercycle 
with pair-connected species, that is. 


Ik 


0 


Ik + Ik+1 - > Ik + 2/fc+i 


(25) 


with fc = l,...,7V — 1, plus closing locm reactions In — > 
0 and In + h In + 2/i (see Figj^. 

In this case, the extra hard-core repulsive reactions 
that need to be imposed are provided by the reactions: 

f (2/, + Ij,h + 21 j) A /, + Ij 
\ {2Ii -f Ij, li 21 j) — > (2/i, 21 j) fori ^ j ^ k 

[_ /i + Ij + Ik —^ li + Ij 


Upon application of Doi-Peliti formalism, the Ho and 
Hi contributions are obtained: 


N N 

Ho = -D +(5 (i)k4>k 

k^l k^l 

N-1 

Hi = —r (j)k+i{^k+i + i)(0fc + 

fc=l 

+ ^i(^i + l)(‘^Af + ( 26 ) 

note that the (jik fields here are not shifted to their non¬ 
trivial vacuum expected values. 

From the formal linear structure provided by the 
hamiltonian dynamics above, it can be argued that fol¬ 
lowing the previous steps leads to a similar conclusion 
regarding the role of dimensionality for the N-member 
hypercycles. However, it is crucial to point out that, even 
though we can probe the system using a coarse-graining 




















scheme as before, both initial conditions and interme¬ 
diate configurations impose a major constraint for this 
kind of systems. This can be understood via two differ¬ 
ent mechanisms. On the one hand, it is well known that, 
for TV > 4, deterministic dynamics of coexistence in hy¬ 
percycles is given by oscillatory behaviour [37H39j . These 
oscillations drive the system close to its extinction state, 
whence, intrinsic-external noise or spatially induced fluc¬ 
tuations may force the system into its inactive phase. 

Even though trying to compute statistical measures 
using (26) is a task beyond the scope of this article, it 
can be conjectured that, provided that sufficiently good 
initial conditions are given, then surface dynamics stands 
as an evolutionary selective advantage for these type of 
hypercyclic systems in contrast with higher dimensional 
layouts. 


V. CONCLUSIONS 

Several studies concerning the origins of life and coop¬ 
eration among molecular replicators have approached the 
problem using a physics-like perspective [101144) . By re¬ 
ducing the system’s complexity to its minimal logic rules, 
we can search for unifying principles of emergent dynam¬ 
ics and universality. 

In order to overcome the limitations imposed by error- 
prone replication in a prebiotic scenario, it has been sug¬ 
gested that a cooperative loop of coupled reactions, the 
hypercycle, is a necessary condition to achieve higher lev¬ 
els of complexity [2 [315]. Similarly, cooperative dynam¬ 
ics might be a requirement to sustain complex ecosys¬ 
tems [13]. The role played by spatial constraints as well 
as stochastic effects has been previously considered by a 
number of authors m, indicating that it can play a key 
role in making possible the stabilisation of hypercycles, 
particularly in relation with their resilience against para¬ 


sites and other sources of disturbance. These studies have 
indicated that in two dimensions the hypercycle can be 
more robust, developing spatial structures that protect 
them (in particular) parasites. It is worth mentioning 
that the study of deterministic hypercycles in three di¬ 
mensions (using reaction-diffusion equations) reveals that 
structures become unstable thus making the resulting dy¬ 
namics less likely to persist [35] . 

In this paper we have analysed a spatially-extended, 
reaction-diffusion field theory of hypercyclic replicators. 
The main target of our study is the simplest, symmetric 
two-member hypercycle involving two cooperative part¬ 
ners. The model includes both second-order reactions 
based on interactions among components of each class as 
well as a linear degradation term. We have studied the 
general conditions for hypercycle survival as provided by 
the RG flow equations. It has been shown that the hyper¬ 
cycle has robust properties (escaping from extinction) for 
dc = 2 dimensions, whereas for d > dc a, dramatic change 
happens, with a much less reduced survival of the molec¬ 
ular cooperators. These results support the view that 
early life required a surface context to persist and evolve 
complexity. Further analysis based on field theoretical 
approaches should consider the role of asymmetries in 
the hypercycle [35| as well as an explicit consideration of 
parasites. Similarly, other classes of replication dynam¬ 
ics (based on templates) as well as additional components 
such as error tails [47] will be explored elsewhere. 
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